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Mass-number and isotope dependence of 
the local microscopic optical potential for polarized proton scattering 
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We derive the local microscopic optical potential U systematically for polarized proton scattering 
at 65 MeV, using the local-potential version of the Melbourne p-matrix folding model. As target 
nuclei, we take ®He and neutron-rich Ne isotopes in addition to stable nuclei of mass number 
A — 4-208 in order to clarify A and isotope dependence of U. The local potential reproduces 
the experimental data systematically and has a geometry similar to the phenomenological optical 
potentials for stable targets. The target density is broadened by the weak-binding nature and/or the 
deformation in unstable nuclei. For the real part of the spin-orbit potential, the density broadening 
weakens the strength and enlarges the radius, whereas for the central potential it enlarges both the 
strength and the radius. The density-broadening effect is conspicuous for halo nuclei such as ''He 
and ^^Ne. Similar discussion is made briefly for proton scattering at 200 MeV. 

PACS numbers: 25.40.Cm, 24.10.Ht, 24.70.+S 



I. INTRODUCTION 

Understanding of nucleon-nucleus (NA) and nucleus- 
nucleus (AA) interactions are one of the most impor- 
tant subjects in nuclear physics. The interactions (opti- 
cal potentials) are necessary to describe not only elastic 
scattering but also more complicated reactions. Actu- 
ally the interactions are inputs of theoretical calculations 
such as the distorted- wave Born approximation (DWBA) 
and the continuum discrctizcd coupled-channels method 
(CDCC) [1H3 to analyze inelastic scattering and transfer 
and breakup reactions. 

In general, NA elastic scattering is less absorptive and 
hence more sensitive to nuclear interior than AA scatter- 
ing. Furthermore, one can determine not only the central 
part but also the spin-orbit part of NA interaction, when 
the incident particle is polarized. From this point of view, 
systematic measurements of polarized proton (p) elastic 
scattering were made so far for stable target nuclei. As 
an outcome of the measurements, some global NA optical 
potentials have been constructed phenomenologically; sec 
for example Refs. [l-H]- In most of the cases, the poten- 
tials are assumed to be local and thereby quite practical 
in many applications such as DWBA and CDCC calcu- 
lations. 

Elastic scattering of unstable nuclei on a polarized pro- 
ton target is a quite good tool of investigating nuclear 
structures of unstable nuclei and interactions between 
a proton and unstable nuclei, but the measurement is 
not easy. The first measurement on the vector analyz- 
ing power Ay was made by Uesaka and his collaborators 
for p+^He elastic scattering at 71 MeV with the inverse 
kinematics |8|. The (/-matrix folding model yields a rea- 
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sonable agreement with the measured vector analyzing 
power [8|. The nuclear- medium effect is thus important 
for the observable. Very lately the effect of ^He breakup 
on the elastic scattering was investigated with the eikonal 
approximation [9| . The medium effect is more significant 
than the breakup effect for Ay. 

The g-matrix folding model is a powerful tool to de- 
scribe NA and AA interactions [l(]| - |l9j . Particularly 
when the Melbourne g-matrix nucleon-nucleon (NN) in- 
teraction [131 is used, the model is successful in reproduc- 
ing polarized proton elastic scattering from stable target 
nuclei systematically with no free parameter |17| . In the 
model, effects of target excitation and breakup in addi- 
tion to the medium effect are taken into account within 
the local-density approximation. The microscopic optical 
potential derived by the model is nonlocal and thereby 
not so practical in many applications. It is, however, 
possible to localize the potential with the Brieva-Rook 
approximation |l2{ . Recently the validity of the approx- 
imation was shown in Ref. [20| . 

From a theoretical viewpoint based on the multiple 
scattering theory [2U]23], a multiple NN scattering se- 
ries in AA collision 1231 is more complicated than that 
in NA collision [2l|, l23]. In this sense, microscopic un- 
derstanding of AA interaction is relatively more difficult 
than that of NA interaction. Very recently, however, for 
AA elastic scattering the local microscopic optical poten- 
tial was derived from the Melbourne g-matrix interaction 
by using both the folding model and the Brieva-Rook lo- 
calization |24| - [2a |. When the local-potential version of 
the Melbourne g-matrix folding model was applied to the 
12q_,_12q scattering at 6.2, 135 and 250 MeV/nucleon 
and ^°-32]\fg_^i2(^ scattering at 240 MeV/nucleon, the 
resultant local potentials reproduced the measured elas- 
tic cross sections and reaction cross sections ctr with no 
free parameter |24| - [26| . In the calculation, the densities 
of unstable nuclei ^''"■^^Ne were evaluated by antisym- 
metrized molecular dynamics (AMD) [23, [2^ with the 
Gogny-DlS interaction [23|, whereas the phenomenolog- 



ical density [301 '^^^ taken for ^^C. The AMD wave func- 
tions suecessfuUy describes low-lying spectra of Ne iso- 
topes [231 ■ The microscopic approach shows that •^°~^^Nc 
in the "island of inversion" have large deformation and 
particularly "^^Ne is a deformed halo nucleus with spin- 
parity 3/2~ [24l - l26l |. The local-potential version of the 
Melbourne g-matrix folding model is thus a powerful 
tool of determining the structure of unstable nuclei. The 
success of the local-potential version of the Melbourne 
g-matrix folding model for proton scattering from sta- 
ble nuclei and '^'■'^^'^Ne+^^C scattering implies that the 
model is reliable also for proton scattering from unstable 
nuclei such as Ne isotopes. 

In this paper, we derive the local microscopic optical 
potential U systematically for polarized proton scatter- 
ing at 65 McV, using the local-potential version of the 
Melbourne ^-matrix folding model. This local nature of 
U makes it easy to clarify global properties of U for both 
stable and unstable targets. As the targets, we consider 
^He and neutron-rich Ne isotopes in addition to stable 
nuclei of mass number A = 4-208. Throughout the global 
analyses, we clarify A and isotope dependence of U. The 
target density is broadened much in ^He and neutron-rich 
Ne isotopes compared with stable nuclei with the same A 
by the weak-binding nature and/or the nuclear deforma- 
tion. We clarify how the density broadening affects the 
shape and the strength of U. This analysis is interesting 
particularly for ^He and ^^Ne, since ^He is a typical two- 
neutron halo nucleus and ^^Ne is a one-neutron halo nu- 
cleus with large deformation. Similar analyses are briefly 
made for proton scattering at 200 MeV. 

We recapitulate the local-potential version of the Mel- 
bourne p-matrix folding model in Sec. Inland show numer- 
ical results in Sec. lIIII Section ITVI is devoted to summary. 



II. THEORETICAL FRAMEWORK 

Proton clastic scattering can be described by the one- 
body Schrodingcr equation 



(Th + ?/-£;)*(+) = 0, 



(1) 



with the proton optical potential U, where E is the in- 
cident energy of proton and T^ is the kinetic energy op- 
erator concerning the relative coordinate R between an 
incident proton and a target. The optical potential U 
can be divided into the central (CE), the spin-orbit (LS) 
and the Coulomb (Coul) component: 



U ^UcE + Ui^sL -(T + Vc 



oul- 



(2) 



In the (7-matrix folding model, U consists of the direct 
and exchange parts |31| . The exchange part is nonlocal, 
but it can be localized with the Brieva-Rook approxima- 
tion Il2i . The validity of the approximation is shown in 
Ref. [20[ . We take this approximation in the present pa- 
per, since the local nature of U makes it possible to make 
global analyses of U over stable and unstable targets. In 



addition to these merits, the local microscopic optical po- 
tential is quite useful in many applications, particularly 
when a potential between stable and unstable nuclei is 
needed. 

The central part Uce is then described by the sum of 
the direct component UP^ and the localized exchange 
component U^^ [l2l. ITsl. Il8l| : 



Uce = I^CE + iWcE = U^k + U, 
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where Vce (Wce) is the real (imaginary) part of Uce and 
s = r — i? for coordinate r of an interacting nucleon from 
the center-of-mass (cm.) of target. The local momentum 
hK{R) = yj2^.j^{E - Uce - Vboui) present in Eq. ([5]) is 
obtained self-consistently, where fij^ is the reduced mass 
of the proton-target system. 

Usually the direct part g^^ and the exchange part g^^ 
of the ^'-matrix interaction are assumed to be a function 
of the local density pa = Pa{r — s/2) at the midpoint 
of the interacting nucleon pair. The direct and exchange 
parts are described as 
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in terms of the spin-isospin component g^^ of the g- 
matrix interaction. 

Similar derivation is possible for the spin-orbit part 



Uj. = K 
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with 



C/™(i?)=-^E jR-sp^ir)9Eipo.i-s;P^)dr, (9) 
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where Vls (W^ls) is the real (imaginary) part of Uls 
and see Appendix [K\ for the definition of p!["'^(i?, s) and 

6q {R,s). As the g- matrix interaction, we take the 
Melbourne interaction [I^l constructed from the Bonn- 
B nucleon- nucleon potential |32| . 

As the target density p, we consider three kinds of 
matter densities. For all the targets, the matter densi- 
ties are evaluated by spherical Hartrec-Fock (HF) cal- 
culations with the Gogny-DlS interaction [23 in which 



the spurious cm. motion is removed in the standard 
manner [2^. For lighter stable nuclei of ^ < 40, we also 
use the phenomcnological proton-density [3^ determined 
from the electron scattering in which the finite-size effect 
of the proton charge is unfolded in the usual way |33|. 
The neutron density is assumed to have the same geome- 
try as the corresponding proton one, since the deviation 
of the neutron root-mcan-squarc (RMS) radius from the 
proton one is about 1% in the spherical HF calculation. 

For ^°~'^^Ne, the matter densities are evaluated also 
by AMD [13, ill with the Gogny-DlS interaction ^. 
This provides the deformed matter densities. The AMD 
wave functions successfully describes low-lying spectra 
of Ne isotopes [13 and ctr of ^^-^^Ne-h^^C scattering at 
240 MeV/nucleon (24|- [26| . and it is free from the spurious 
cm. motion; see Ref. [2^, [23] for the details of AMD 
calculation. 

Particularly for '^^Ne with quite small one-nucleon sep- 
aration energy, we have to do a sophisticated calculation 
to make a tail correction to the AMD density, since the 
AMD density is inaccurate in its tail region. In principle, 
the ground state of ■^-'^Ne can be described by the ■^"Ne-f n 
cluster model with core ('^"Ne) excitations, and one can 
solve coupled equations for the model with the resonating 
group method (RGM) in which the ground and excited 
states of "^°Ne were constructed by AMD. This was done 
in Ref. [li,[ll|. This method is referred to as AMD- RGM 
in this paper. 

The matter density pimim'if) is obtainable from the 
ground state wave function ^j™^ with spin /, its z- 
component m and parity n as 

Pi.ni^'{r) = ($^:ri E'^(^^ - ^ - ^)i*g"'")' (11) 

i 
21 

^Y.PHir)iIm'Xp\Im)Y:^if), (12) 



A=0 



where X is the cm. coordinate and r^ is the coordinate 
of i-th nucleon, and the summation of A in Eq. (J12p runs 
for even numbers. As a standard manner in the folding 
model, only pjj~ is taken. This procedure is exact for 
even nuclei with / = 0, but approximate for odd nuclei 
with / 7^ 0. It is, however, confirmed in Ref. [2^ that 
the procedure is a quite good approximation. 



III. RESULTS 

In this section wc analyze polarized proton scattering 
a.t E = 65 and 200 McV, using the local-potential version 
of the Melbourne _g-matrix folding model, although the 
analysis is mainly focused on the case oi E — 65 MeV. 
In the folding model, we also consider the Franey-Love 
i-matrix NN interaction [SJ] . As target nuclei of proton 
scattering, we consider stable nuclei of ^He, ^^C, ^^O, 
20Ne, 24Mg, 40Ca, ^^Zr, ^ospb and unstable nuclei of '^Hc 
and neutron-rich Ne isotopes. 



A. Proton scattering from stable nuclei at 65 and 
200 MeV 

In this subsection, we mainly consider proton scatter- 
ing at 65 MeV and discuss the case of 200 MeV at the 
end of this subsection. 

Figure [1] shows the elastic cross section da/<Kl and the 
vector analyzing power Ay as a function of the scattering 
angle ^c.m. for proton scattering bX E = 65 MeV from 
^He, i^C, 40Ca and ^ospb targets. The solid lines denote 
results of the local microscopic potential derived from the 
Melbourne g-matrix interaction. The results are consis- 
tent with those [17[ of the nonlocal microscopic poten- 
tial derived from the Melbourne g-matrix interaction, al- 
though the target densities used are different between the 
two calculations. The g-matrix folding model (solid line) 
yields better agreement with the experimental data |35l - 
|42| than the i-matrix folding model (dashed line). The 
medium effect shown by the difference between the solid 
and dashed lines are more significant for Ay. For lighter 
targets of ^He, ^^C and ''"Ca, the agreement of the g- 
matrix folding model with the data is better for the phe- 
nomcnological density (solid line) than for the HF density 
(dotted line). However, the g-matrix folding model with 
the HF density still keeps reasonable agreement with the 
data. The HF density based on the Gogny-DlS inter- 
action |29| is thus quite useful to derive the microscopic 
optical potential U systematically. 

Figure [5] shows reaction cross sections a-n as a func- 
tion of A. Closed circles (cross symbols) denote results 
of the g-matrix (t-matrix) folding model with the phe- 
nomcnological target density for A < AQ and the HF 
target density for A > 40. The g-matrix folding model 
yields better agreement with the data than the i-matrix 
folding model. The results of the g-matrix folding model 
agree with the data [43l - |47l | within 10% error in a wide 
range of 4 < ^ < 208. Since a^^ is sensitive to Wce, the 
part is thus derived with at least 90% accuracy. 

As a merit of deriving the local microscopic optical 
potential [/, wc can directly compare U with the phe- 
nomcnological optical potentials. Here wc take two types 
of the phenomcnological optical potentials. One is the 
Koning-Delaroche (KD) optical potential [j| based on 
the standard potential search in the framework of the 
Schrodingcr equation, and another is the Dirac phe- 
nomenology optical potential [5|, [6| based on the poten- 
tial search in the framework of the Dirac equation. As 
an example, we compare the three optical potentials for 
the p+4"Ca system at S = 65 MeV in Fig. H The 
microscopic optical potential is calculated from the phe- 
nomenological target density. The imaginary spin-orbit 
part Wls is about ten times as small as VlSj so that 
it hardly affects the present elastic scattering. We then 
compare the three optical potentials for Vqe, Woe and 
Vhs- The KD optical potential (dashed line) has sim- 
ilar R dependence to the Dirac phenomenology optical 
potential (dotted line) . In the peripheral region, further- 
more, R dependence of the two phenomcnological optical 
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FIG. 1: (Color online) Angular distributions of the elastic 
cross section do /d£l (a) and the vector analyzing power Ay (b) 
for proton scattering at _E = 65 MeV from ''He, ^^C, *°Ca and 
■^'''^Pb targets. The solid (dotted) lines represent the results of 
the g-matrix folding model with the phenomenological (HF) 
density, whereas the dashed lines correspond to results of the 
t-matrix folding model with the phenomenological density. 
Experimental data are taken from Refs. [35| - |43 |. 
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FIG. 2: (Color online) A dependence of reaction cross sec- 
tions (7r at iJ ~ 65 MeV. Closed circles (cross symbols) stand 
for results of the g-matrix (t-matrix) folding model with the 
phenomenological density for A < 40 and the HF one for 
A > 40. Experimental data are taken from Refs. [43h47| 1. 



where E =47.9 MeV for *He, 65.5 MeV for ^^C, ^®0, ''°Ca, 
2°*Pb, 47.0 MeV for ^^Ne, 
for 'J^Zr. 



48.0 MeV for ^"Mg and 60.8 MeV 



several percent. 
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FIG. 3: (Color online) R dependence of the microscopic opti- 
cal potential (solid line), the Koning-Delaroche (KD) optical 
potential (dashed line) [J] and the Dirac phenomenology op- 
tical potential (dotted line) 0, [y] for p+'^'Ca scattering at 
£ = 65 MeV. 



potentials is similar to that of the microscopic one (solid 
line). At small i?, the deviation of the microscopic op- 
tical potential from the phenomenological ones is small 
for VcE, but becomes sizable for Vls and Wce- For Fls 
the deviation little affects the measured da/dfl and Ay 
up to 90°, but for Wce the deviation slightly enhances 
(Tr and consequently overestimates the measured ctr by 



For the p+'^He system, the g-matrix folding model docs 
not reproduce the measured clastic cross section at back- 
ward angles 9c. m. ^ 90°. A possible origin of this devi- 
ation is the heavy-ion (■^H) exchange process that is not 
included in the folding model. The process generates a 
parity-dependent potential [43 . We then rewrite Vce as 



VcE^VcE(l + a(-l)^) 



(13) 



to incorporate the effect in the folding model phcnomeno- 
logically. The parameter a is fitted to the experimental 
data [35|, |3^, and the resulting value is a = 0.13. In 
Fig- SI the solid and dashed lines show da/dVt and Ay cal- 
culated by the (7-matrix folding model with a = 0.13 and 
0, respectively. The effect of the parity-dependent poten- 
tial is appreciable at intermediate angles 60° < Oc.m. ^ 
90° and significant at backward angles at 0c. m. ^ 90°. 
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FIG. 4: (Color online) Effect of the parity-dependent poten- 
tial on the elastic cross section da/dQ (a) and the vector an- 
alyzing power Ay (b) for p-l-^He scattering at i? = 65 MeV. 
The solid (dashed) line stands for the result of the g-matrix 
folding model with a = 0.13 (0), where the phenomenological 
density is taken for *He. Experimental data are taken from 
Refs. m,^- 



For the p+^^C system, the measured Ay is not repro- 
duced around 6'c.m. = 100° by the g-matrix folding model. 
Even the phenomenological optical potentials [^-Q can- 
not reproduce it. The parity-dependent potential does 
not solve this problem, since the potential enhances the 
elastic cross section largely at Oc.m. > 90°. Therefore this 
is an open problem. Except for the measured Ay around 
^c.m. = 100° for the p+^^C system, the g-matrix folding 
model yields good agreement with the data for targets of 
A = 4-208. 

Now we briefly analyze proton scattering at _E := 
200 MeV as an example of intermediate energies and 
a typical case of beam energies in Radioactive Isotope 
Beam Factory (RIBF). 

Figure [5] shows angular distributions oida/dfl and Ay 
for the scattering from stable nuclei "^He, ^^C, ""^Ca and 
208pj-, rp]_^^ definition of lines is the same as in Fig. [1] 
The g-matrix folding model with the phenomenolo gica l 
density, shown by the solid line, reproduces the data |5ll - 



|54| particularly at forward angles 6'c.m. ^ 60°. (p+^^C 
65MeV : 9 = 128°,200MeV : 9 = 60°)When Fig. 
[5] is compared with Fig. [1] as a function of the trans- 
ferred momentum, 0c. m. = 60° in Fig. [5] corresponds to 
6'c.m. = 128° in Fig. [f| Therefore, 6'c.m. > 60° in Fig. [5] 
correspond to very backward angles 0c. m. > 128° in Fig. 
[T] In this sense, it is natural that in Fig. [S]the solid lines 
are deviated from the data at 0c. m. ^ 60° particularly for 
lighter targets of ^ < 12. Even at 0c. m. ^ 60° for lighter 
targets of ^ < 12, the g-matrix folding model with the 
phenomenological density (solid line) yields better agree- 
ment with the data than that with the HE density (dot- 
ted line). The medium effect, shown by the difference 
between the solid and dashed lines, is still significant at 
E = 200 MeV. 



B. p+^'He scattering at 71 and 200 MeV 

As a typical example of proton scattering from unsta- 
ble targets, we first consider p-t-^He scattering at 71 MeV, 
since the experimental data is available not only for 
da/dO. but also for Ay. The ^He density is obtained by 
the ^He-t-n-l-n orthogonality-condition model; see Refs. 
[49I |50 | for the detail. As the ^He-density in the model, 
we take either the phenomenological or the HE density. 
In Eig.ini the solid (dashed) line shows da/dfl and Ay cal- 
culated by the (^-matrix (f-matrix) folding model with the 
phenomenological ''^He density, whereas the dotted line 
corresponds to the result of the g-matrix folding model 
with the HE '''He density. The 5-matrix folding model 
(solid line) yields better agreement with the data @, l55| 
than the t-matrix folding model (dashed line) . The solid 
line is consistent with a result of the nonlocal microscopic 
optical potential derived from the Melbourne g-matrix 
interaction [8| , although the target densities used are dif- 
ferent between the two calculations. 

In Fig. [6l the difference between the dashed and dotted 
lines comes from that between the phenomenological and 
HE ''He densities. Comparing Fig. IH] with Fig.[TJ one can 
see that the difference of *He density yields smaller effect 
on da/dfl and Ay in the p-|-^He scattering than in the 
p-|-^He scattering. The difference of ^He density is thus 
partly masked by the density of extra two neutrons in 
the p+^He scattering. 

In Fig. [6l the nuclear medium effect is shown by the 
difference between the solid and dashed lines, and the 



effect is such more significant for Ay. The difference 
mainly comes from the fact that Wce calculated with 
the g-matrix folding model is much less absorptive than 
that with the t-matrix folding model; see Fig. [7] for R 
dependence of U. 

As mentioned above, the g-matrix folding model yields 
better agreement with the data than the t-matrix folding 
model. The agreement is, however, not perfect particu- 
larly for Ay. The parity-dependent potential little works 
at forward angles < 90° in which the experimental data 
is available. We then make the optical potential search 
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FIG. 5: (Color online) Angular distributions of the elastic 
cross section da jdQ. (a) and the vector analyzing power Ay 
(b) for proton scattering eX E = 200 MeV from *He, ^^C, 
''"Ca and ^°**Pb targets. See Fig. [T]for the definition of lines. 
Experimental data are taken from Refs. [5ll - [54l [. 



to see what causes the difference between the theoretical 
and experimental results. Here we change the strength 
of U only as 






Coul 



(14) 



with adjustable parameters {j\f J\^ J\f J)^) ■ The 
best fit is obtained when {!^^ JW ^ JW" ^ Jw) == 
(0.7,0.7,0.9,4.0), but dajdO. and A^^ at 6* > 40° are not 



FIG. 6: (Color online) Angular distributions of the elastic 
cross section da /dO, (a) and the vector analyzing power Ay 
(b) for p+'^He scattering at £; = 71 MeV. See Fig. [l] for the 
definition of lines. Experimental data are taken from Ref. [^, 

in. 



sensitive to changes of f^ and fy^ and hence these pa- 
rameters are not determined sharply. 

In Fig. [51 the result (solid line) of the potential search 
is compared with the original result (dashed line). Small 
reduction of da/dfl from the dashed to the solid line at 
6 > 50° mainly comes from reduction of Vce, and large 
change of Ay from the dashed to the solid line at 9 > 40° 
is mainly originated in large enhancement of IVls • Thus 
the large enhancement of Wls is necessary, but it is not 
clear what causes the large enhancement of VFls ■ This is 
an interesting question to be solved in future. 

Now we briefly analyze proton scattering from ^He at 
E = 200 MeV. Experimental data is not available for this 
energy at the present stage, but the measurement of Ay 
is planned in RIBF. Figure [9] shows angular distributions 
of da/dO. and Ay for proton scattering from ^He at i? = 
200 MeV. The definition of line is the same as in Fig. [Bl 
As shown in Fig. [SJ the results of the g-matrix folding 
model (solid lines) are consistent with the data for stable 
target nuclei. Hence the model prediction, shown by the 
solid line in Fig. [HI is expected to be reliable also for ^He 
particularly at forward angles 0c. m. ^ 60°. At backward 
angles 0c. m. ^ 60°, the solid line is different from the 
dashed and dotted lines. This means that inclusion of 
the medium effect and usage of the phenomenological 
''He density seem to be necessary there. 
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FIG. 7: (Color online) R dependence of microscopic optical 
potential for p+^'He scattering at E = 71 MeV. Solid (dot- 
ted) line denotes the potential of the g-matrix folding model 
with the phenomenological (HF) ^He density. Dash line is the 
potential of the t-matrix folding model with the phenomeno- 
logical *He density. 
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FIG. 9: (Color online) Angular distributions of the elastic 
cross section da/dQ (a) and the vector analyzing power Ay 
(b) for proton scattering from "^He at £ = 200 MeV. See Fig. 
|6]for the definition of lines. 
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FIG. 8: (Color online) Potential search for p-f^He scat- 
tering at E = 71 MeV. The calculations are done with 
the p-matrix folding model with the phenomenological *He 
density. The solid line represents the best-fit result with 
ifv^Jw^Jv^Jw) = (0.7,0.7,0.9,4.0), whereas the dashed 
line corresponds to the original result. Experimental data are 
taken from Ref . M, Im . 



C. A dependence of microscopic optical potential 

The f/-inatrix folding model is successful in reproduc- 
ing the data on polarized proton scattering for targets 
oi A = 12-208. In general, the (/-matrix folding model 
is considered to be less accurate for light targets such 
as ^'^He, since the g-matrix interaction is evaluated in 
nuclear matter and the local density approximation used 
seems to be more reliable for heavier targets. Neverthe- 
less, the agreement of the model calculations vifith the 
data is still reasonable for light targets ''■^He. We then 
derive global properties of the local microscopic optical 
potential for proton scattering at 65 McV. The potential 
is calculated with the g matrix folding model from the 
phenomenological target density for A < AO and the HF 
density for A > 40. 

Since the microscopic optical potential docs not have 

any simple shape, we consider the volume integral J and 

1 

the RMS radius {R'^) '^ for each part of U: 



J ^ F(R)dR, (R^)'- =J R^F{R)dR/J, (15) 

where F{R) represents each part of U. As shown in Figs. 
|3]and[7l Wls is much smaller than Vls and hence yields 
no significant effect on the present systems. Furthermore, 
the meaning of the volume integral and the RMS radius 
is not clear for VFls, since Wls is an oscillating function 
of R. For these reasons, we do not consider the volume 
integral and the RMS radius of Wls here. 



Figure [10] shows the RMS radius of L'^ as a function 
of A. For each part of U, A dependence of the RMS 
radius is fitted for stable targets oi A = 4-208 by a curve 
aiA + a2A^/^ + a3A^/^ + a4; see TableUfor the parameter 
set (01,02,03,04). The fitting fine is referred to as the 
stable-nucleus line. For each panel, the lower figure shows 
the ratio of each RMS radius to the stable-nucleus line. 
Since ^He has a halo structure, the ^He RMS radius is 
enlarged by the structure from the stable-nucleus line by 
~18% for VcE, ~25% for Wqe and -12% for V^s- Large 
enhancement of the ^He RMS radius in Wqe means that 
the reaction cross section is a good tool of searching light 
halo nuclei and determining their radii. 
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FIG. 10: (Color online) A dependence of the RMS radius of 
the microscopic optical potential at 65 MeV. Three panels cor- 
respond to VcE, WcE and Vhs- See Table [I] for the parameter 
set of the stable-nucleus (dotted) line. 





ai 


02 


03 


04 


VcE 


-0.0023 


0.0104 


0.8508 


1.0578 


WcE 


-0.0076 


0.0473 


0.8732 


1.332 


Vls 


-0.0033 


0.0137 


1.2305 


-0.0682 



TABLE I: Parameter fitting for A dependence of the RMS 
radius of U. All the parameters are presented in units of fm. 



We consider the reason why the enhancement of ^He 
RMS radius is larger in Vce than in Vls- For this pur- 
pose, we show A dependence of RMS radii of p and dp/ dr 
in Fig. 111! In the short-range limit of the g-matrix in- 
teraction, the shape of Vce (^ls) is approximately deter- 
mined by that of p {dp/dr). Enhancement of ^He RMS 
radius from the stable-nucleus line is smaller for dp/dr 
than for p, as shown in Fig. [TT] see Table HI] for the 
parameter set of the stable- nucleus line (dotted line) . 
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FIG. 11: (Color online) A dependence of the RMS raidus of 
p (upper panel) and dp/dr (lower panel). See Table HIl for the 
parameter set of the stable-nucleus (dotted) line. 



Figure [T2| shows the volume integral of t/ as a function 
of A. The volume integral is also enhanced at A = 6 
from the stable nucleus line by the halo structure of ^He. 
This enhancement is larger in Wqe than in Vce- In Vls, 
meanwhile, the volume integral is suppressed at A = 6 
from the stable nucleus line. The halo structure of ^He 
makes p broader and hence suppresses dp/dr, and conse- 
quently weakens the strength of Vls- 





ai 


02 


a-i 


aA 


p 

dp/dr 


-0.0074 
-0.0079 


0.0378 
0.0402 


0.9578 
1.2478 


-0.0022 
-0.2322 



TABLE II: Parameter fitting for A dependence of tlie RMS 
radii of p and dp/dr. All the parameters are presented in 
units of fm. 





ai 


02 


as 


a4 


VcE 
WCE 

Vls 


301.28 
81.082 
0.093 


-59.695 

92.201 

6.24 


272.22 
301.22 
95.645 


-125.78 
-439.22 
105.11 



TABLE III: Parameter fitting for A dependence of the volume 
integral of U . All the parameters are presented in units of 
MeV fm^ 



Finally we show A dependence of (Tr in Fig. [T3] for 
proton scattering at £' = 65 MeV. The reaction cross 
section is enhanced for ^He by ^ 35% from the stable- 
nucleus line (dotted line); see Table ITVl for the parameter 
set. This enhancement is a result of the corresponding 
enhancement in the RMS radius and the volume integral 
of WcB- 

Throughout all the analyses in this subsection, one can 
find that the broadening of ^He density due to the halo 
structure (weak-binding nature) yields significant effects 
on U . The effect is different between the central and 
spin-orbit forces. For the central force, it enhances the 
strength and the RMS radius by about 20% ~ 40%. For 
Vls, meanwhile, the broadening enlarges the RMS radius 
by about 10% but weakens the strength by about 20%. 



D. Isotope dependence of microscopic optical 
potential 

The density broadening due to the weak-binding na- 
ture appears not only in ®Hc but also in the general un- 
stable nuclei where the one-nucleon or the two-nucleon 
separation energy is small. In addition to the weak- 
binding nature, the nuclear deformation also makes the 
diffuseness of density larger, i.e., induces the density 
broadening. 

The success of the local-potential version of the Mel- 
bourne g-matrix folding model for ^°~^^Ne-|-^^C scatter- 
ing and proton scattering from stable nuclei at 65 MeV 
may indicate that the model is reliable also for p-|-^°~'^^Ne 





ai 


a2 


a-j, 


a4 


O-R 


-0.6797 


45.534 


122.34 


-170.78 



TABLE IV: Parameter fitting for A dependence of ctr. All 
the parameters are presented in units of mb. 
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FIG. 12: (Color online) A dependence of the volume integral 
of the microscopic optical potential at 65 MeV. See Table IIIII 
for the parameter set of the stable-nucleus (dotted) line. 



scattering at 65 VIeV. We then analyze this scattering 
with the model in order to investigate the weak-binding 
and the deformation effect on U and ctr. The analysis is 
interesting particularly for ■^^Ne, since it is a one-neutron 
halo nucleus with large deformation. 

Figure [T3] shows isotope dependence of a^i for 
p_l_20-32-[\^g scattering at 65 MeV. The dotted line denotes 
the stable- nuclei line for ctr. The cross symbols represent 
results of the ij-matrix folding model with the spherical 
HF densities, whereas the closed circles correspond to 
those with the AMD densities in which Ne isotopes are 
deformed. The deviation of the cross symbols from the 
dotted line represents the weak-binding effect, whereas 
that of the closed circles from the cross symbols shows 
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FIG. 13: (Color online) A dependence of the reaction cross 
section for proton scattering at i5 = 65 MeV. See Table IIVI 
for the parameter set of the stable-nucleus (dotted) line. 



FIG. 14: (Color online) Isotope dependence of the reaction 
cross section for p+ ^ Ne scattering at 65 MeV. See Table 
IIVI for the parameter set of the stable-nucleus (dotted) line. 



the deformation effect. These effects are more apparent 
in the lower panel where ctr for Ne is normalized by 
that for stable nuclei with the same A. At 20 < ^ ^ 23, 
the deformation effect is more significant than the weak- 
binding effect. Since the weak-binding effect becomes 
large as A increases, the effect is more important than 
the deformation effect at 24 < A < 30. For A = 31,32, 
Ne isotopes have no bound state in the spherical HF cal- 
culation. In this sense, both the deformation and the 
weak-binding effect are important there. 

For ^"'^Ne, the one-neutron separation energy is quite 
small. We then do a RGM-AMD calculation to make a 
tail correction to the AMD density. The closed square 
denotes a result of the RGM-AMD density, and the de- 
viation of the closed square from the stable-nucleus line 
represents the net effect of the halo structure and the 
deformation. This effect is detectable from the measure- 
ments of CTR or the one-neutron removal cross section. 
Actually, this was done for the "^^Nc + ^^C system at 240 
MeV/nucleon [H-ii, [H, [13 • 

Figures [TS] and [12] show the volume integral and the 
RMS radius of U for p+^^'^^Nc scattering at 65 MeV, 
respectively. The weak-binding nature and the nuclear 
deformation play the same role in U as in ctr. Actually, 
the deformation effect is more important than the weak- 
binding effect at 20 < A < 23, whereas the latter is 
more significant than the former at 24 < A < 30. The 
density broadening due to the weak-binding nature and 
the deformation in Ne isotopes yields the same effect on U 
as that due to the halo structure (weak-binding nature) in 
^He. For '^^Nc as a typical case, it suppresses the volume 
integral by about 10% and enlarges the RMS radius by 
about 5% for Vls, whereas it enlarges both the volume 
integral and the RMS radius by about 10% for Vce and 
M^CE- The imaginary central part Wqe is more sensitive 
to the density broadening than the other parts. Thus ctr 
is a good quantity to investigate the density broadening. 



IV. SUMMARY 

We have derived the local microscopic optical potential 
systematically for polarized proton scattering at 65 MeV, 
using the local-potential version of the Melbourne g- 
matrix folding model that was successful in reproducing 
nucleus- nucleus scattering [24 265 ■ The local microscopic 
optical potential for proton scattering well reproduces the 
experimental data on stable targets systematically and 
has a geometry similar to the phenomenological optical 
potentials such the Koning-Delaroche [j| and the Dirac 
phenomenology 0, Q one. Also for p-|-^He scattering at 
71 MeV, the g-matrix folding model yields better agree- 
ment than the i-matrix folding model. Therefore we can 
say that the local-potential version of the Melbourne g- 
matrix folding model is a reliable model to describe both 
nucleon-nucleus and nucleus-nucleus scattering. 

The local microscopic optical potential U is also quite 
useful in many applications. Furthermore, the local na- 
ture makes it easy to clarify global properties of U over 
both stable and unstable target nuclei. As target nuclei, 
we have then considered ^He and neutron-rich Ne iso- 
topes in addition to stable nuclei of mass number A = A- 
208. Particularly, it is interesting to compare properties 
of U in ^Hc and '^-'^Ne with general properties of U in sta- 
ble nuclei, since ^He is a typical two-neutron halo nucleus 
and '^-'^Ne is a one-neutron halo nucleus with large defor- 
mation. Throughout the global analyses, we have clari- 
fies A- and isotope dependence of U. The target density 
is broadened in ^Hc by the weak-binding nature and in 
neutron-rich Ne isotopes by both the weak-binding na- 
ture and the nuclear deformation. For the central part of 
U the density broadening enlarges both the strength and 
the root-mean-square radius, but for the real spin-orbit 
part of U it enlarges the root-mean-square radius but 
weakens the strength. The broadening effect is conspicu- 
ous for halo nuclei such as ^He and ^^Ne. The statement 
mentioned above is true also for intermediate energies 
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FIG. 15: (Color online) A dependence of the volume integral FIG. 16: (Color online) A dependence of the RMS radius 



of the microscopic optical potential for p+^''~'^^Ne scattering 
at 65 MeV. See Table IIIII for the parameter set of the stable- 
nucleus (dotted) line. 



of the microscopic optical potential for p+^''~'^^Ne scattering 
at 65 MeV. See Table [I] for the parameter set of the stable- 
nucleus (dotted) line. 



such as 200 MeV. 
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Appendix A: Definitions of quantities 



The explicit forms of Sq {R, s) and Pi {R, s) are 
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